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Small particles advected in a fluid can collide (and therefore aggregate) due to the stretching 
or shearing of fluid elements. This effect is usually discussed in terms of a theory due to Saffman 
and Turner (J. Fluid Mech., 1, 16-30, (1956)). We show that in complex or random flows the 
Saffman- Turner theory for the collision rate describes only an initial transient (which we evaluate 
exactly). We obtain precise expressions for the steady-state collision rate for flows with small Kubo 
number, including the influence of fractal clustering on the collision rate for compressible flows. For 
incompressible turbulent flows, where the Kubo number is of order unity, the Saffman- Turner theory 
is an upper bound. 
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1. Introduction. Suspensions of small particles in a 
fluid, such as aerosols or colloids, are ubiquitous in the 
natural world and in technology. Such systems may be 
unstable due to collisions of the suspended particles giv- 
ing rise to aggregation or to chemical reaction. These 
collision processes are fundamental to understanding the 
formation of raindrops from clouds [H, or the growth of 
planets from dust in a protostellar environment 0] . 

Collisions between particles will always occur due to 
Brownian diffusion. The rate of collision may be greatly 
increased if the fluid is in macroscopic motion Q , such as 
turbulence [J]. If the inertia of the suspended particles 
is significant, the suspended particles may move relative 
to the fluid and this effect can greatly increase the rate 
of collisions @, @, 0] . However if the particles are suffi- 
ciently small that the effect of their inertia is negligible 
(advective limit), the particles may be brought into con- 
tact by the effect of shearing motion in the flow. This 
process was discussed by Saffmann and Turner [J] and 
their theory has been the starting point for most subse- 
quent work on advective collisions in smooth flows. Here 
we point out a previously unremarked assumption in [J] , 
which implies that their expression for the collision rate 
describes a short-lived transient. In the case of incom- 
pressible flow their expression is an upper bound: we give 
the first precise results on the steady-state collision rate 
of their model. Both the Saffman- Turner theory and our 
expression are compared with numerical simulations in 
Fig. [T£i,b (for synthetic incompressible flows in two and 
three dimensions respectively; the details are explained 
later). In both cases we observe an initial transient given 
by the Saffman- Turner theory, and at larger times the 
collision rate approaches a significantly lower value which 
is in agreement with our theory. 

Reference 0] treats the local flow in the vicinity of a 
particle as if it were a steady hyperbolic flow, whereas 



in reality it must fluctuate. We argue that the appro- 
priate theory for the collision rate must depend upon a 
dimensionless parameter termed the Kubo number, Ku. 

2. Summary. In this letter we give the first precise 
expressions for advective collisions in smooth flows. We 
obtain asymptotic results in the case of flows with small 
Kubo number, in both two and three dimensions, in- 
cluding flows with Richardson diffusion and compressible 
flows where particles cluster onto a fractal set. We also 
describe new results on the evaluation of the Saffman- 
Turner expression for the initial collision rate in two and 
three dimensions. 

3. A formula for the collision rate. We start by dis- 
cussing a general formula for the collision rate, which is 
very similar to that proposed by Saffmann and Turner 0] . 
We model processes where particles coalesce upon con- 
tact with probability unity. The particles collide when 
their separation falls below 2a (where a is the particle 
radius) . For simplicity we initially discuss two spatial di- 
mensions. Let v r (r,9,t) be the radial component of the 
velocity of the fluid relative to a test particle at (r, 0) 
(in polar coordinates). We make the approximation that 
the velocity field is not perturbed by the suspended par- 
ticles: this is valid in the Epstein low-pressure regime 
Q and is a reasonable approximation Q in conventional 
fluids. The collision rate for a single particle is then the 
rate at which points distributed randomly with density 
no cross into the disc of radius 2a centered on the origin. 
The instantaneous collision rate may then be written as 

R=-2an d6v r (2a,6,t)Q(-v r {2a,6,t))x(2a,6,t) (1) 
Jo 

where the Heaviside function 0(-) selects for particles 
entering a disc of radius 2a, no is the number density of 
the suspended particles in the vicinity of the test par- 
ticle, and xi r i 0, t) is zero if the fluid at (r, 9) at time t 
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FIG. 1: Collision rate as a function of time for particles ad- 
vected in an random flow, a incompressible flow in two spa- 
tial dimensions, Uo = 1, f] = 0.1, no — 10 3 , r = 0.001 (o) 
and r = 0.004 (□), a = 0.003. Theory for initial transient 
(|11|) . dashed lines, for collision rate ©, full lines. Error bars 
correspond to one standard deviation, b three spatial dimen- 
sions, uo = 1, r = 0.004, n = 0.1, a = 0.02, and n = 2xl0 3 . 
Theory for initial transient (|18|l , for collision rate (|17|l . c two- 
dimensional compressible flow. Parameters: Uo = 1, rj = 0.1, 
r = 0.004, a = 0.001, n = 10 3 , and V = 2. Initial transient 
(|16|) . dashed line; colllision rate (|15[) , solid line, d steady hy- 
perbolic flow in the vicinity of a test particle, e shows region 
(in rest frame of test particle) occupied by particles which will 
collide with test particle after time t — 20r (parameters as in 
Fig. la, □), f same but after time t = 10 3 r. 

has previously passed through the disc of radius 2a, and 
unity otherwise. The collision rate of a single particle R a 
is the time-average of ([IJ, yielding an exact and general 
expression for the advective collision rate, but in the gen- 
eral case it is difficult to evaluate since the characteristic 
function \ depends on the history of the flow. 

The relative speed v r is determined as follows 0]: for 
small particles a<?;we approximate the flow u(r, 0) in 
the vicinity of the test particle by a linearised equation 
for the small distance X from the test particle: X = 
A(0, 0)X, where A(r, t) is the strain matrix of the flow 
u(r, t). We then have v r = X ■ n((9) where n(#) is a unit 
vector normal to the surface at a point with polar angle 
6. At the point X = 2an(9) on the surface of the disc 

v r = 2an • An . (2) 

Eqs (|1I2D appear as the key expressions for the collision 
rate in but without the factor \- Ref- describes 



the flow in the vicinity of a test particle as if it were a 
steady hyperbolic flow (Fig. [TJi). In that case, particles 
advected by the local flow only pass through the circle 
r = 2a once (Fig. so that the factor \ has no effect. 
In general, however, the flow in the vicinity of the test 
particle will fluctuate, as illustrated in Fig. QJ, so that the 
directions of the eigenvectors of A rotate. The function 
%(r, 9, t) equals zero for regions where the flow at r = 2a 
has reversed direction, so that points that were originally 
advected away return to the surface. This reduces the 
collision rate below the Saffman- Turner estimate. 

There is one limit in which the Saffman- Turner the- 
ory is exact: this is the limit of very short times, when 
the initial condition ensures that x(r, 9,0) = 1. For later 
times, it might be expected that the Saffman- Turner the- 
ory is an upper bound on the collision rate. This is true 
if the flow field is incompressible. For compressible flows 
there is a complication, because particles may cluster to- 
gether, and the density in the vicinity of a test particle 
may be larger than the average density, no- 

Although it is not possible to make a precise theory for 
the steady-state collision rate in the general case, this can 
be done in the limit where the velocity field of the fluid 
is rapidly fluctuating. The effect of fluctuations depends 
on a dimensionless parameter termed the Kubo number. 
A random or turbulent flow may be characterised by a 
correlation length rj, a correlation time r, and a charac- 
teristic velocity uq. The Kubo number Ku = uor/r) is 
a dimensionless measure of how rapidly the flow fluctu- 
ates. It might be argued that the Saffman- Turner theory, 
which treats the flow as if it were steady, would apply to 
flows with very large Kubo number. However, such flows 
arc unphysical, because the velocity field transports spa- 
tial fluctuations so that they become temporal fluctua- 
tions, implying that the Kubo number is never large. In 
the case of fully-developed turbulence, Ku — 0(1). 

4. Collision rate for small Kubo number. In the limit 
Ku — * (that is, where the correlation time approaches 
zero), the relative separation of two particles undergoes 
a diffusion process determining the collision rate. 

We consider the case of two spatial dimensions. By ho- 
mogeneity we can assume that the particles are at and 
X = (Xi, X2) and by isotropy we can take X2 = and 
write X = (r, 0). The fluctuation SXi of the separation 
in a short time 6t satisfies (6Xi) = 0. The correlation of 
these increments is (SXi6Xj) = 2Dij5t, with 

Dij = \ f dt([ Ul (X,t)- Mj (O,t)][ Uj (X,0)- Uj (O,0)]). 

" J -00 

(3) 

These 'diffusion constants' are computed as follows. 
We consider a velocity field u(r,t) generated by a 
'stream function' ip(r,t): u — V A ipiiz, where n z is 
a unit vector in the z-direction, with correlation function 
(ip(r,t)ip(0,0)) = C(r,t), where r = \r\. For X = (r,0) 
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we find Di 2 = D 2 \ — and 



dt [C"(0,t)--C'(r, t)\ 



/oc 
dt [C"(0,t)-C"(r,t)} 
-OO 



(4) 



where the prime denotes a derivative w.r.t. r. 

We set up a Fokker- Planck equation for the separation 
r of two particles. The change in the value of r in St 
is Sr = y[x + Sx) 2 + Sy 2 — x = Sx + Sy 2 /(2x) , so that 
(Sr) = D 22 (r)St/r and (Sr 2 ) = (Sx 2 ) = 2D lx (r)St. Wc 
obtain the Fokker-Planck equation 



dp d 
dt dr 



D 22 (r)p d , v 

+ -^-(D n (r)p) 

r or 



(5) 



The collision rate is determined by a steady-state solution 
of ([5]) with a constant flux — J. The boundary condition 
is p(2a) = 0, which corresponds to particles which reach 
r = 2a being absorbed, so that there are no re-collisions. 
The collision rate is i? a = J. To enforce the boundary 
condition at r = 2a, it is sufficient to consider an approx- 
imate solution valid at small values of r. Expanding the 
diffusion constants for small r we find Du(r) ~ Vr 2 and 
D 22 (r) ~ 3£>r 2 , with 
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V = - dtC""(0,t) 



we obtain the desired steady-state solution 
p(r) = J/(2Vr) + Ar . 



(6) 
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To determine the constant A we match the solution (JT]) 
for J = to the exact steady-state zero-flux solution of 
([5]). The latter satisfies 



d log p 1 

dr r 
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(8) 



From ([4]) we find a(r) = 0. The density is thus pro- 
portional to r, normalisation gives p(r) — 2imor. Wc 
conclude that A = 2-kuq in (0). The collision rate is 
determined by setting p(2a) = in ([7|) and solving for 

Ra.= J: 



R a = 16ttD n a 2 



(9) 



5. Richardson diffusion. In a steady, fully developed 
turbulent flow particle separations exhibit Richard- 
son diffusion d(r 2 )/dt oc r 4 / 3 in the inertial range r n 
fioj ] (where 77 is the Kolmogorov length). This implies 
Dij(r) oc r 4 / 3 . It is of interest to consider the collision 
rate in a flow exhibiting Richardson diffusion in the limit 
Ku — > 0, to establish whether additional complications 
arise when considering a multiscale flow (although real 
turbulence has Ku = O(l), and is at the borderline of 



applicability of our theory). In the following we argue 
that (HI) is unchanged for a flow exhibiting Richardson 
diffusion, such that Dij(r) ~ r* in the inertial range. 

At first sight, given eq. ([5|), this power-law scaling of 
the diffusion constants suggests that there should be a 
power-law behaviour of p(r) in the inertial range, with 
an exponent which depends upon £. This would be in- 
consistent with the requirement that p(r) ~ 2irnor for 
r 3> rj (as it must be in a homogenous and isotropic 
flow). However, given equations (jlj, we see that C(r,t) 
contains a term which can be approximated by a multi- 
ple of r 2+ £ for r in the inertial range. It follows that the 
coefficients of D\\ and D 22 are related, D\\(r) — Cr^ 
and D 22 (r) = C(l + £) r ^- Comparison with ([5|) shows 
that p ~ r in the inertial range, as expected. Since the 
velocity field in two dimensions is generated by a stream 
function, it follows immediately that a(r) — 0, and thus 
p(r) — 2-Kn^r for r <C rj as well as for r ^> 77. We con- 
clude that our result j9]) for the collision rate of particles 
advected in incompressible flows at small values of Ku is 
insensitive to what happens in the inertial range. 

6. Computer simulations. The numerical simulations 
shown in Figs. [1^ and b employ a Gaussian random 
stream function with spatially and temporally station- 
ary, spatially isotropic statistics with mean value zero 
and C(r,t) = (u\r\ 2 j2) exp[-|f|/r - r 2 /(2n 2 )}, for which 
V = «§t/(2?7 2 ). The particles are initially randomly po- 
sitioned, and they are regarded as having collided when 
their separation falls below 2a. Our plots show the rate 
of collision with a test particle which is also advected 
by the velocity field. Upon collision with the test par- 
ticle, colliding particles are removed. The collision rate 
was determined by dividing by t the number of particles 
which collided in the time interval [0,t]. This yields ac- 
curate results for the initial transient, Rq, and for the 
steady-state collision rate, R a . Only a small fraction of 
particles had collided by the end of each run (less than a 
few percent). Fig. []Ji shows the collision rate as a func- 
tion of time for particles advected in a two-dimensional 
incompressible flow for small Kubo number. We observe 
good agreement with the analytical result ©, full line. 
In the following we show how to derive the initial tran- 
sients shown in Fig. [T^, dashed lines. 

7. Exact result for the initial transient. The initial 
transient is obtained by averaging equations ([T|) , ([2]) over 
the ensemble of random strain matrices determined by 
the statistics of the stream function ^(0,0). In order 
to compute the average of R, we decompose the matrix 
A which appears in ^ into a symmetric Aq and an 
antisymmetric Aj] part and sort the eigenvalues <Ji,i — 
1, . . . , d of the symmetric part in order: o\ < o 2 <,..., < 
<j c i. In two dimensions we obtain for a given A 



R(cr 2 ) = 8noa 2 a 2 



(10) 



We average over the distribution of the largest eigenvalue 
of Aq, P(<t 2 ) = {2a 2 n 2 /ul)exp[-(rji] 2 /ul], where a 2 > 
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0, and find for the initial rate of collision, Rq: 
Rq = i^/nnoa 2 uq/t] . 



(11) 



This initial transient is shown in Fig. [T]a (dashed lines). 
Comparing with ([5]), we have i? a ~ KvlRq, demonstrating 
that the effect of eliminating multiple collisions is signifi- 
cant when the velocity field u(r,t) is rapidly fluctuating. 

8. Collision rate in the presence of fractal clustering. 
When the velocity field has a compressible component, 
the particles cluster onto a fractal set This may al- 



ter the collision rate, because particles which cluster to- 
gether are expected to collide more frequently [f|. Here 
we present the first exact results on this effect. When 
considering compressibility effects we restrict ourselves 
to two spatial dimensions, because the natural experi- 
mental realisation of particles advected by compressible 
flows involves particles floating on the surface of a turbu- 
lcnt fluid, such as the experiments discussed in 
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We consider a partially compressible flow of the form 
u = (V A tph z + pV(j>)/y/l+/P where if) and 4> are 
independent Gaussian random functions with correlation 
function C(r,t) and (3 > 0. As before we seek a solution 
of ([5]) with constant flux — J. We have 



D 2 2{r)-~ 



-oo 
oo 



dt[C"(0,t)- 



r(l+/3 2 ) 1+/3 2 
C"(r,t) f3 2 C'{r,t) 
1+P 2 r(l+/3 2 ) 1 



so that A22/A.1 T = ((3 2 + 3)/(3/3 2 + 1) as r 
The steady-state solution of ([5]) is for r <C 77 



p{r) = [(T+1)/(T-1)] J/{AVr)+Ar 



r-2 



0. 



(13) 



The constant A is determined, as before, by matching to 
the small-r behaviour of the exact zero-flux solution of 
([5]) with the diffusion constants (fT"2"l) : 



p(r) = 2irnor exp 



(14) 



dr'a(r') 



In order to match (fT3|) with (fT4|) we evaluate (|14[) for 
small values of r. One obtains p(r) ~ 2irnoi](r/ri) D2 ~ 1 
where D2 — T — 1 is the correlation dimension of the 
fractal set onto which the particles cluster (this dimen- 
sion was first determined in Ref. [HI]). Matching with 
(fl~3f we obtain (for 1 < T < 3) 



R a = [(r - l)/(r + 1)] 8ttV n oV 2 (2a/ V ) 1 



(15) 



For the initial transient in compressible flows we find: 

Rq = (8V^/Vl + r) n a 2 u /n . (16) 

Comparison with (fl5|) shows that eq. ([I]) no longer gives 
an upper bound of the collision rate upon setting x — 
1. This is because Eqs. Q an d (fill) assume a uniform 
distribution of particles. But for times t ^> |Ai + A2 1 — 1 



(here Ai > A2 are the Lyapunov exponents of the flow), 
the particles cluster on a fractal set [ll[ which can lead 
to an increase in the collision rate (Fig. lc). 

There is a critical point at T = 1 where Ai becomes 
negative and path-coalescence occurs [3] ■ Experimental 
evidence indicates that V = 1 for surface flows above 
a turbulent fluid At this critical point we obtain 

R a = 4irVn n 2 \og(2a/n). 

9. Three spatial dimensions. We consider a three- 
dimensional incompressible flow generated by three ran- 



dom vector fields [15( . The calculation of R a proceeds in 



a fashion similar to the derivation of 



R a = 967rX>n a J 



V = 



ld 2 D n W 



dr 2 



(17) 



As in the two-dimensional case we find that Richardson 
diffusion does not alter the result p7|) . The initial tran- 
sient in three dimensions is derived as outlined above. 
We have obtained an explicit expression for R(<J2, 03) in 
terms of the largest eigenvalues a% , 03 of Aq , but could 
not perform the average of 02, 03 analytically. Numerical 
averaging for our Gaussian-correlated model gives 



Rq w 24 noa uo/rj . 



(18) 



Results ([T7]) and (fTS)) are shown in Fig. QJd in comparison 
with numerical simulations. 
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